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Abstract 

We consider numerically the growth of neutrino asymmetry in active-sterile neutrino os- 
cillations in the early universe. It is shown that the final sign of the asymmetry can be 
highly sensitive to small variations of the oscillation parameters. We find regions which are 
completely or partially chaotic, but also regions where the sign remains very robust. The 
consequences for atmospheric neutrino oscillations and primordial nucleosynthesis are then 
discussed. In the completely chaotic region the predicted 4 He-abundance has an inherent 
arbitrariness AY ~ 10 -2 . 



Active-sterile neutrino oscillations in the early Universe is a fascinating possibility with 
far-reaching consequences e.g. for nucleosynthesis [1-12] and CMB radiation [|H|. Nucle- 
osynthesis considerations in particular have made it possible to place stringent constraints 
on model building aimed at understanding the observed neutrino anomalies in terrestrial 
observations. In the very first papers |l|] it was observed that the mixing with an active 
species (S77(2)-doublet) endows the sterile (S77(2)-singlet) neutrino with effective interac- 
tions, which can be strong enough to bring the sterile species in equilibrium. The ensuing 
excess energy density would result in a failure of the nucleosynthesis explanation of the ob- 



served light element abundances (14| . This line of reasoning was put to a solid computational 
foundation in refs. [Q, |J, and these results were later reproduced in ref. ||. Of particular in- 
terest was the observation that nucleosynthesis is in conflict with v^ — //^-oscillation solution 
to the atmospheric neutrino problem ]7]]. 

Already in f|| it was noted that nucleosynthesis constraints [§], || depend on the reasonable 
assumption that the leptonic asymmetries are not many orders of magnitudes larger than 
the baryonic one. More specifically, e.g. for a mass squared difference \§m 2 \ = 10 -4 eV 2 a 
large initial asymmetry, L™ > 1CT 5 (here L a = (N a — Na,)/N 7 ) would suppress the effective 
mixing angle so much that the equilibration would never take place. This observation was 
later revived by Foot and Volkas 0, who basing on this and another previously observed 
effect, an exponential growth of leptonic asymmetry ||, suggested an interesting way to 
circumvent the nucleosynthesis constraints without invoking unnatural initial conditions [|10j . 
Their scenario assumes a novel mass-mixing scenario, where a v T — zvrnixing, with carefully 
chosen parameters, produces a large leptonic asymmetry (but does not equilibrate u s ), which 
suppresses the subsequent v^ — ivmixing angle and thereby prevents the ^-equilibration from 
taking place. Some details concerning the growth of the asymmetry in this scenario are still 
under debate 0, 112" . 



It was later observed by Shi (jTa] that the period of exponential growth exhibits chaotic 
features and therefore, while the amplitude of the final asymmetry is robust, its sign appeared 
to be essentially arbitrary. This raises some interesting questions: for example, is the sign of 
L u sensitive to the fluctuations in the initial conditions, like in the baryon asymmetry? If so, 
one should expect a large suppression in the effective asymmetry present at the important 



epoch for the z/ M — z/ s -oscillations due to diffusion effects. For this purpose it is important 
to establish the extent of possible chaotic or regular regions in the parameter space. In this 
letter we have studied the dependence of the sign of L v on neutrino mixing parameters. 
We find a rather clearcut division of the parameter space into non chaotic and partly or 
completely chaotica regions. In chaotic regions the final sign of the asymmetry is indeed 
found to be highly sensitive also to fluctuations in the initial conditions. 

Another, more direct consequence follows from the fact that sign(L) affects the computed 
4 He abundance, either directly in the case of v e — v s oscillations, or when induced by large 
L v created in v^ — v s or v T — v s oscillations and later transferred to z/ e -sector via active-active 



oscillations ||16[ . It then follows that possible chaotic behavior will constrain our chances to 
draw any definite conclusions about the effects of sterile neutrinos on Big Bang nucleosyn- 
thesis, as we will discuss below. 

In the early Universe neutrinos experience frequent scatterings, which tend to bring their 
distributions into thermal equilibrium. The requisite mathematical formalism is therefore 
very different from the one particle approach valid for description of accelerator physics 
(beams) and even solar neutrinos. Indeed, the objects of interest are the (reduced) density 
matrices for the neutrino and antineutrino ensembles 

p„ = ip (l + P), p p ^Ip (l + P). (1) 



is 



Solving full momentum dependent kinetic equations for p v (p) and pp(p) 0, [17|, [18], |19 
obviously a very difficult task. Instead, we employ the momentum averaged equations for 
P = P((p)), with (p) ~ 3.15T, which should be expected to give a good approximation 
for the full system H. (Our preliminary studies with full momentum dependent equations 
support this assumption). Moreover, for the parameters we are interested in, one can neglect 
the collision terms so that Pq remains a constant and can be set to a unity. The coupled 
equations of motion then are (for definiteness we shall focus here on v T — v s oscillations; 



1 We do not claim here that the system exhibits chaoticity in the mathematical sense of the definition of 
chaos; we merely mean that the system is sensitive to small variations of the parameters. 



other cases are obtained from this by simple redefinition; 

P = V x P - DP T 

P = v x P - DP T (2) 

where P = dP/dt and we defined P^ = -P^x + P y y. In the case of v T — v s oscillations the 
damping coefficient is D ~ 1.8G F T 5 |||, and D ~ D to a very high accuracy. It is convenient 
to decompose the rotation vector V as 

V = V x ±+(V + V L )z, (3) 

with the components 

V x = A sin 26 

V = -Acos26 + 5V T 

V L = V2G F N^L, (4) 

where 6 is the vacuum mixing angle, A = 5m 2 /2(p) and the photon number density iV 7 = 
2((3)T 3 /tt 2 . The effective asymmetry L appearing in the leading contribution Vl to the 
neutrino effective potential is 

L=~L n + L v . + L Vlt + 2L Vr {P) (5) 

where n refers to neutrons and we have assumed electrical neutrality of the plasma. The 
remaining piece to the effective potential SV T is given by 0, |2T 



5V T = -V2G F N y A T { ^, (6) 

with A T = 14£(4)/£(3) ~ 12.61. The rotation vector for antineutrinos is simply V(L) = 
V(— L). The coupling of particle and antiparticle sectors occurs through the asymmetry 
term, where 

L Vr {P) = l£ + l{P M -P,). (7) 

with L™ being the initial v T asymmetry. 

2 Interested reader can find these redefinitions for example from B. 



Even the simplified one state-quantum kinetic equations (||]) are very difficult to handle 
numerically, because of the vast difference in the time scales involved (Hubble expansion 
rate, matter oscillation frequency and the width of the resonance, for example) on one hand 
and due to extremely strong coupling induced by the asymmetry term on the other. The so- 
called static approximation employed in || reduces the system to one first order differential 
equation. Unfortunately it is not really suitable for the treatment of oscillations at the 
resonance, since many of its basic assumptions - that the system is adiabatic, that the 
MSW-effect can be neglected, and that the rate of change of lepton number is dominated by 
the collisions - break down at the resonance. 

These considerations emphasize the need for a very careful numerical approach. In prac- 
tice, the accuracy is much improved if one makes separation between the large (~ number 
density) and small components (~ asymmetry) in equations (Q). To this end we change the 
variables into 

Pt = Pi±Pi, (8) 

in terms of which (0) become 

P+ = -V P y + -V L P y -DP,+ 

P+ = V P+ + V L P--DP+ 

Pt = V X P+ 

P- = -V P y --V L P y + -DP x T 

P" = V P- +V L P: - DP~ 

P; = V X P- (9) 



We have studied numerically the behaviour of the system described by (|D|) as a function of 
the oscillation parameters 5m 2 and sin 2 29, and in particular the evolution of the asymmetry 
(0). Of crucial importance in this evolution is the occurrence of the resonance at Vq = if 
Sm 2 < 0. Inserting the appropriate parameters, one finds that the resonance temperature is 
given by g 

T rcs ~ 16.0 (\5m 2 \ cos26) 1/6 MeV, (10) 

where Sm 2 is given in units eV 2 . Far above the resonance the damping terms tend to 



suppress the off-diagonal elements Py and moreover, the system is driven towards the 
initially stable fixed point L = 0. As soon as the system passes the resonance however, 
L = becomes an unstable fixed point and two new locally stable and degenerate minima 
corresponding to the solutions oi Vl + V = appear; these are given by the condition 
\L\ ~ ll(|5m 2 |/eV 2 ) (MeV/T) 4 . 

The system is roughly analogous to a ball rolling down a valley that branches to two, 
passing via a saddle configuration. Initially the branching into these two new valleys can be 
very shallow and it may stay that for a long time. Once on the side of the bifurcation, the 
ball still keeps on passing over the central barrier (the continuation of the old stable fixed 
point to an unstable extremum) until the barrier grows too high, or until friction (damping) 
reduces energy enough, and it gets trapped to one of the valleys. It is easy to picture in 
ones mind that in a case of a very shallow bifurcation a small change in the initial conditions 
(L m ), or in the shape of the valleys (oscillation parameters), can very much affect which 
minimum the system finally chooses to settle in. 

In Fig. 1 we show the asymmetry L Ut fl7|) resulting from solving equations (H) for two 
sets of parameters. In Fig. la the resonance is rather narrow and only one oscillation occurs 
before the system is trapped into the minimum with a negative sign of L Ut . The subsequent 
oscillations about this new local minimum are quickly washed away by the damping terms. 
In contrast, Fig. lb shows an example of oscillation parameters with which the bifurcation 
into new local minimum is extremely slow and for a long time there is hardly any barrier 
between the two minima with opposite signs of L Vt , and the system oscillates thousands of 
times before settling down to a minimum with positive L Vt . After settling down, the further 
evolution of the asymmetry follows a power-like behaviour. These results agree well with 
those of ref . [|TJ] . 

It is instructive to look a little more carefully into how the system approaches the res- 
onance. Before the resonance the off-diagonal P ± components are very near zero and P~ 
near the value L = 0. Just before the resonance the P+ components begin to increase, 
which triggers both the growth of P~ components and the decrease of P 2 + . As Vq changes 
the sign at the resonance it creates an instability in the equation for P~, which eventually 
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Figure 1: The evolution of the z/ T -asymmetry in the neigbourhood of the resonance. In Fig. 
la 5m 2 = -KT 1 and sin 2 26 = 2 x 10" 9 ; in Fig. lb 5m 2 = -KT 2 and sin 2 26 = 2 x 10" 5 . 

strongly pushes P~ to negative direction. The simple coupling of P~ to P~ in (f|) then 
drags P~ along leading to a rapid growth of L Vt . So far these phenomena have not much 
affected the evolution of P + variables (which have continued to grow). Eventually however, 
the exponential growth of P~ terms causes the Vl term in the equations for P^ and Pj~, 
insofar neglible, become dominant. Large Vl then forces P% and Py to change sign and grow 
to opposite direction until Vl again changes the sign. Additionally, the ensuing oscillatory 
motion of the P+ components induces the oscillation into other variables as well, leading to 
the exponentially large oscillation pattern observed in L Vt . 



To find out the extent of chaotic and/or regular behaviour of sign(L), we have scanned 
through the parameter space depicted in Fig. 2, which shows the sign of the final asymmetry 
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Figure 2: The distribution of the final sign of the neutrino asymmetry in the mixing pa- 
rameter space. Negative sign(L) is plotted in black, positive sign(L) in white. The initial 
asymmetry was chosen to be L m = 1CT 10 . 



L with the initial value L m = 10~ 10 . As can be seen, the structure of sign(L) is highly 
complex. In the upper left hand corner, extending downwards to large 9, there is a regular 
region with no change in the asymmetry. Its existence is relatively easy to understand: 
this is the region where only the very first oscillation is carried out before the sign of the 
asymmetry is fixed. Since the direction of the first oscillation is determined by the sign of the 
initial asymmetry L m (not necessarily the initial iv-asymmetry) , the sign of final neutrino 
asymmetry in this part of the parameter space should indeed be regular and fully determined. 
The bands seen in the left hand side of Fig. 2 are formed as the system goes through two 
or more oscillations. In this region the number of oscillations is slowly increased as 9 grows 
leading to less determined sign(L) but it still can hardly be described as chaotic yet. 

In addition to the two more or less regular regions there are regions where sign(L) appears 



to be chaotic. The interval 10~ 2 < \5m 2 \ < 1 contains a very complicated structure. For 
\5m 2 \ > 1 one may discern some tendency for positive L to prevail, while the region with 
\5m 2 \ < 10 -2 appears to be pure white noise. In the lower right-hand corner of Fig. 2 
(below the gray line), with 1 5m 2 1/ sin 2 2$ < 10, oscillations in L will continue past the 
neutrino freeze-out and will not settle into any definite value. The boundary of the regular 
region above which L is positive is given approximately by 

/ -5m 2 \ 2 ( sin 2 2#\ 2 , N 

log ? -0.44 log — ~-l. 11 

V & io- 5 - 8 ev 2 ; V 10 / 

We conjecture that the chaotic behaviour occurs only when the oscillating period is 
long. We have also explored in finer detail a restricted region of parameters in the chaotic 
regime, without finding any structure. It is possible however, that the final sign of L in 
this region is affected by the accumulated numerical error originating from the extremely 
high number of oscillation periods. In this sense proving a true chaoticity is of course not 
possible. Nevertheless, if the system is sensitive to numerical error, it should be expected 
to be sensitive to parameter fluctuations as well, so the general pattern of rapid sign(L) 
fluctuations is expected to be robust. Finally, the region in the upper right-hand corner does 
not correspond to a large asymmetry, but it is merely the region where v s is fully equilibrated 
||, |7|, |j and the absolute value of final L is very small. 

Changing the initial value for L m does not change the picture qualitatively, although the 
structures evident in Fig. 2 shift slightly to the left when L m is increased. Moreover, the 
changes saturate at L in > 10 -9 . These changes, or their absence, are very interesting how- 
ever: In Fig. 3 we have plotted the value of L at the temperature T = T res — 2.5, as a function 
of L m for three representative choices of parameters. The first set, with sin 2 29 = 10~ 7 and 
5m 2 = — 10 4 eV 2 corresponds to the stable region with positive L in Fig. 2. As expected, L 
remains positive independently of the initial value L m . In fact the dependence turned out to 
be smooth and linear showing that not only is the sign robust, but also that the numerical 
solution is well under control. The second set, with sin 2 29 = 10~ 6 and Sm 2 = —1 eV 2 , lies 
in the intermediate region where positive L predominates, and the same dominance is seen 
as a function of the inital value L m . The last set with sin 2 29 = 10~ 6 and 5m 2 = — 10~ 3 
eV 2 corresponds to the chaotic region. It is evident that sign(L) is very sensitive to initial 
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Figure 3: The sign of the asymmetry L Vt at T = T rcs — 2.5 MeV as a function of the initial 
asymmetry L m for three sets of parameters (sin 2 26, 5m 2 /eV 2 ): a) (10~ 7 ,— 10 4 ), b) (1CT 6 ,— 1) 
and c) (l(r 6 ,-l(r 3 ). 

conditions, displaying clear randomness as a function of L m . 



The final value of sign(L) has consequences for both atmospheric neutrinos and primor- 
dial nucleosynthesis. It has been proposed that Super-Kamiokande results for atmospheric 
neutrinos, which lie in the forbidden zone [§], |||, 0, might still allow a active-sterile mixing 



solution if the asymmetry growth is taken into account |TTj . Although the oscillation param- 
eters in the case of atmospheric neutrinos are in the region where asymmetry growth is not 
expected, it has been argued that other neutrino oscillations could induce a large asymmetry 
in the active-sterile sector which the oscillations cannot damp. 

If the outcome of neutrino oscillations is highly chaotic, the validity of such a scenario 



9 



might be suspect. However, we found a large region in the parameter space where sign(L) is 
very robust with respect to small variations of the mixing parameters. No chaoticity should 
be expected there with respect to other small perturbations, such as local perturbations in 
L m , either. It is in these stable domains where one would expect that the mechanism of ref. 
|T0| can be successful. 

In the region where sign(L) is chaotic in the oscillation parameter space, it was also found 
to be sensitive to fluctuations in L m ; these are predicted to be generated for example during 
the QCD phase transition, or in scenarios of electroweak baryogenesis [22, [23), 24]. In such 



case causally disconnected regions would be expected to develop large asymmetries with a 
random sign distribution. It has been argued that then the nucleosynthesis constraint on 
active-sterile mixing would be even more stringent, because of additional MSW conversion 
taking place in the boundaries of domains with different sign(L) |23J. However, our results 
indicate that the new constraints obtained in |2J| may be overly optimistic, because for a large 
part of their excluded region we have found sign(L) to be stable against small fluctuations; 
hence in no domain formation should be expected to occur in the first place. 

Determining the sign of L is important also for considering the effect of the electron 



(anti) neutrino spectrum distortions on the light element abundances [jig, ^TJ. When the 
momentum spectrum gets distorted from its thermal equilibrium value the neutron to proton 
freezing ratio will change. Direct u e <-► v s oscillations obviously can induce such distortions, 
but also scenarios where large asymmetry is first generated in v^ r <-» v s and then transferred 
to electron neutrino via v e <-► v^ T oscillation, could have considerable effects on the electron 
neutrino spectrum. It turns out that positive sign(L) has the effect of decreasing and negative 
sign(L) of increasing 4 He abundance [ IB[ , so that the difference is AY ~ 10~ 2 , with some 



dependence on the oscillation parameters. Because the oscillation parameters cannot be 
measured with an arbitrary accuracy, it follows that in the region where the sign(L) is 
chaotic, the role of resonant active-sterile neutrino mixing in Big-Bang Nucleosynthesis can 
not be reliably estimated. Rather, in this region, depicted in Fig. 2, there always remains an 
arbitrariness in the 4 He abundance given by AY ~ 10 -2 , which should be considered as a 
source of systematic error. 

In the region where the sign is stable, more concrete conclusions can be drawn. However, 



10 



in this region L is positive and only a rather small negative shift in the helium abundance 
—0.005 < AF < was found [^6[ for these parameters. Interestingly enough, such a shift 
could ameliorate the apparent conflict of the nucleosynthesis theory viz-a-viz observations 



14 



Our results in this paper are based on an averaged momentum description of the neutrino 
ensemble. Some effects, like the diffusion of the asymmetry between different momentum 
states, would seem to indicate the need for using full momentum-dependent kinetic equa- 
tions. This is rather hard, since one has to deal with exponential growth in every momentum 
state and the width of the resonance is, for most of the parameter space, so small that one 
needs a very large number of bins to complete the task. Our preliminary results with mo- 
mentum dependent kinetic equations support the results presented here. 

This work has been supported by the Academy of Finland under the contract 101-35224. 
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